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Abstract. A pair trade is a portfolio consisting of a long position in one asset 
and a short position in another, and it is a widely applied investment strategy 
in the financial industry. Recently, Ekstrom, Lindberg and Tysk studied the 
problem of optimally closing a pair trading strategy when the difference of the 
two assets is modelled by an Ornstcin-Uhlcnbcck process. In this paper we 
study the same problem, but the model is generalized to also include jumps. 
More precisely we assume that the above difference is an Ornstein-Uhlenbeck 
type process, driven by a Levy process of finite activity. We prove a verification 
theorem and analyze a numerical method for the associated free boundary 
problem. We prove rigorous error estimates, which are used to draw some 
conclusions from numerical simulations. 



1. Introduction 

A portfolio which consists of a positive position in one asset, and a negative 
position in another is called a pair trade. Pairs trading was developed at Morgan 
Stanley in the late 1980's, and today it is one of the most common investment 
strategies in the financial industry. The idea behind pairs trading is quite intuitive: 
the investor finds two assets, for which the prices have moved together historically. 
When the price spread widens, the investor takes a short position in the outper- 
forming asset, and a long position in the underperforming one with the hope that 
the spread will converge again, generating a profit. A main advantage of pairs 
trading is that the short position can, in principle, remove any exposure to market 
risk. For a historical evaluation of pairs trading we refer to [6]. 

To model the pair spread the authors in [3] proposed a mean reverting Gaussian 
Markov chain which they considered to be observed in Gaussian noise. Recently, in 
[2] the authors suggested the continuous time analogue, the so called mean reverting 
Ornstein-Uhlenbeck process. In this paper we generalize the model of the spread to 
also include possible jumps. Let (fi, J 7 , P) be a complete probability space where 
the following processes are defined in such a way that they are independent: 

- A standard Brownian motion W = {Wt}t>o- 

- A Possion process N x = {Nf} t >o with intensity A > 0. 
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- A sequence of independent random variables {X'£}'j*L 1 with common con- 
tinuous symmetric density ip. Moreover, the support of ip is contained in 
the interval (— J, J) for some J > 0. 

Define the compound Poisson process C X:V = {C^' v }t>Q in the usual way as 

k=l 

and denote the filtration generated by W, C >{f and the null sets of _F by F = 
{•^t}t>o- It is well known that this filtration satisfies the usual hypotheses (see 
for example |10j). From now on, when we say that a process is a martingale, 
submartingalc or supcrmartingale we mean that this is with respect to F. 

Let the difference U — {Ut}t>o between the assets be the unique solution of the 
stochastic differential equation 

(1.1) dU t = -fxU t dt + adWt + dC^' p , t > 0, 

where /.i > 0, a > 0. (The solution of equation (|1.1[) is usually called a generalized 
Ornstcin-Uhlcnbeck process or an Ornstein-Uhlenbeck type process.) Sometimes 
we will denote the driving Levy process in (fLTJ) by Z cr ' A:V , i.e. 

ZS> X > V = aW t + C^, t>0. 

As discussed in [2], there is a large risk associated with a pair trading strategy. 
Indeed, if the market spread ceases to be mean reverting, the investor is exposed to 
substantial risk. Therefore, in practice the investor typically chooses in advance a 
stop- loss level a < 0, which corresponds to the level of loss above which the investor 
will close the pair trade. For a given stop-loss level a < define 

(1.2) r a =inf{t>0: U t < a}, 

the first hitting time of the region (-co, a], and the so called value function 

(1.3) V(x) = supE x [U TaAr ] xeR, 

T 

where the suprcmum is taken over all stopping times with respect to U. (Here and 
in the sequel means expected value when Uq = x.) The major interest here is to 
characterize V, and perhaps more importantly, to describe the stopping time where 
the supremum is attained. Since the drift has the opposite sign as U, we have no 
reason to liquidate our position as long as U is negative. On the other hand, if U 
is positive, then the drift is working against the investor and for large values of U 
the size of the drift should overcome the possible benefits from random variations. 
Moreover, since the jumps are assumed to be symmetric, this indicates that there 
is a stopping barrier b > with the property that we should keep our position when 
Ut < b and liquidate as soon as Ut > b. We note that we cannot be sure to close the 
pair trade at any of the boundaries a or &, because the spread can exhibit jumps. 
This was not the case in [2] and it is the major reason for the additional difficulties 
encountered in the present paper. 

General optimal stopping theory (described for example in [H Ch. 3]) leads us 
to believe that the value function is given by V = u, where (u, b) is the solution of 
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the free boundary problem 

Guu{x) = 0, x € (a, 6), 

(1.4) u(sc) = x, x £ (a, 6), 

u'(b) = 1. 

Here £/[/ is the infinitesimal generator of U, which is defined on the space of twice 
continuously diffcrentiable functions / : M —> US with compact support: 

(1.5) guf(x) = ?-f"(x)-vxf'(x) + \[ (f(x + y)-f(x)) l p(y)dy, xtR. 

Moreover, the stopping time where the supremum in (|1.3|) is attained should be 

(1.6) n = inf{i > : E/ t > b}. 
Indeed, our first result is a so called verification theorem. 

Theorem 1.1. Assume that (u,b) is a classical solution of (|1.4[) with 

a) Gu u (x) < 0, for x > b, 

b) u(x) > /or 

TTien m(x) = V(a:) = ~E x [U Ta/ \ Tb ], for where V is given by ()1.3j) . 

Remark: As seen from the assumptions on y>, we are assuming that the absolute 
value of the jumps of the process {Ut}>o are bounded. The reason is that on the 
financial market, an asset cannot jump to arbitrarily large levels. If nothing else, 
the jumps are bounded by all the money in the world. 

The rest of the paper is organized as follows. In Section[2]we prove Theorem ll.il 
and in Section[3]we discuss a numerical solution of the free boundary problem (|1.4j) . 
We also present strong evidence for the existence and uniqueness of a solution of 

(HU). 



2. Proof of Theorem 11.11 

Before we start to prove Theorem 11.11 we need to recall some facts. From the 
general theory in [5] we get that the boundary value problem 

Guu(x) =0, i£ (a, b), 
u{x) — x, x $ (a, b), 

has a unique classical solution and that such a solution belongs to the space 

C 2 (M \ {a, b}) n C 1 (K\ {a, 6}) n C(R). 

Moreover, the finite left and right limits of u' and u" exist at a and b. Although 
these facts follow from [5], we present in Theorem 13. II a self-contained proof for the 
simpler situation that we consider here. Hence, if (it, b) is a classical solution of 
(|1.4[) . then necessarily 

u E C 2 (R \ {a, 6}) n C 1 (R\ {a}) n C(R) 

with finite left and right limits of u' and u" everywhere. Furthermore, recall a gen- 
eralized version of Ito's formula for convex functions (see for example [10l Ch. 4]): 
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Theorem 2.1 (Meyer-Ito formula). Let X = {Xt}>o be a semimartingale and let 
f be the difference of two convex functions. Then 

f{X t ) = f(X ) + f D-f(X s -) dX s 
Jo+ 



(f(X.)-f(X B -)-D-f{X a -)AX a ) 



0<s<t 
f oo 



L\{X)dy.{y), 



where D f is the left derivative of f , fi is a signed measure which is the second 
generalized derivative of f and {L1(X)} t >o is the local time process of X at a. 

Due to the regularity of u it can be written as a difference of two convex functions 
(see Problem 6.24 in [7J Ch. 3]). Moreover, the second derivative measure fi of u 
can be split into two parts fi = \i c + fid, where the continuous part fi c is given by 
dfir = u" dx and the discrete part fid = 6 a is a point mass at a. Here, u"(x) denotes 
the second derivative of u at x except at the points a and &, where it denotes the 
right second derivative (which we know is finite). By Corollary 1 of the Meyer-Ito 
formula in |10] , we can now write 



1 r°° , , , If* „,„ . .„„ 1 



(2.2) 



^LV(U)dfi(y) = - u"(U^)d[U,U} c s + -LnU)(u'(a+)-u'(a-)) 
y f u"{U s -)ds + h?(U){u'(a+) - u'(a-)), 



where [U, U] c denotes the continuous part of the quadratic variation [U, U]. 
Furthermore, by using (jl.ip and the compensated Poisson random measure 

N z (dt,dy) = N z {dt,dy) - \dtip(y) dy, 

where Nz denotes the jump measure associated with Z a ' X > v , we get 

f D-u(U s -)dU s + V (u(U s )-u(U s -)-D-u(U s -)AU s ) 

J °+ 0<s<t 

i [ U s -D-u(U s -)ds + o- [ D-u(U s -)dW s 
Jo Jo 

u(U s - +y)- u{U s -)) N z (ds, dy) 



(2.3) 



o+ 

t 



+ X {u{U s -+y)-u{U s -))tp{y)dyds. 



Summing up, we now have for t > 
u{U t ) = u(U 



J (y/({/ s _) - iOJ.-D-u(U a -)} ds 



(2-4) +xf f (u(U 3 -+y)-u{Us-))<p(y)dyds 



+ lmU)(u'(a+)-u'(a-))+M t , 
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where 

M t = a f D-u(U s -)dW s +11 {u{U s - + y) - u(U 3 -)) N z (ds,dy). 
Jo Jo+ Jr 

Since u is Lipschitz, has a bounded left derivative and since the jumps density iphas 

a finite swe get that {M t }t>o is a martingale. 

Lemma 2.2. Assume a£M and Uq > a. Then a.s. At (U) = for all t > 0. 

Proof. Fix a € K and assume Uq > a. Since the local time process {L°} t >o is 
continuous in t it is enough to prove that for fixed t > we have L% At(^0 = a - s - 
From [H p. 217], we get that 

\L a TaM (U) = (U TaM -ay- hu.->a}(U a - a)" 

0<S<T a At 

fT a At 



/ l{U s _<a}^a- X] !{(7 3 _<a}(^.s - «) + - 

J0+ n^-^_ w 



0<s<r a At 

Futhcrmore, from the fact that U s > a for all < s < r a A t, we get that U s - > a 
for all < s < r a A t and from the left continuity of U s ~, we can conclude that 
we also have U TaA t- > &■ From that and by splitting the integral and the sum, we 
obtain 

l L r a At(U) = l{U TaM -=a}{U Ta M ~ a)~ + l { U TaAt _= a} (U TaAt ~ a) 



- l{u TaAt -=a}{U Ta At ~ a) + + I l{U s -=a}dU s 

- J2 hu.-=*}(U'-a)~ 



0<s<T a M 

-T a f\t 



i>T a /\l — 

/ l {U,-=a} dU s - ^2 1 {U s -=a}(U s - a) + . 

J0+ n^„^_ a + 



0<s<T a M 

From the observation that if £/ s _ = a for some < s < r a A t, then s is a jump time 
and the jump must be in the up direction, we conclude that the right hand side of 
the last expression is zero and so we are done. □ 

Remark: In a similar way one can show that, if a < Uo < 6, then 

L a TaATbM (U) = and L b TaArbAt (U) = for t > 0. 

Proof of Theorem ] Since u(x) = V{x) = E x [U TaATb ] = x, when x < a, we can 
assume that x > a. Define Y t = u(U TaAt ), t > 0. By using (|2.4[) . Lemma I2T21 the 
expression (|1.5p for the generator of U, and Q1.4p . we get 

rr a At 

Y t = u{x) - I At^ s -l{ir s _>6} ds 

+ X / (it(C/ s _ +y) -u(U s -))tp(y)l{ Ua _> b ydyds + M TaAt . 

JO JR 

Property a) and the martingale property of {M TaAt } give that {Y t }t>o is a super- 
martingale. Furthermore, from property b) we get that Y t > U TaA t, for t > 0, and 
since 

(2.6) U TaM > a - J, t>0, 



(i 
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we can apply the optional sampling theorem (see [7]) and obtain 

Ex[U TaAT } < Ex[Y T ] < E X [Y ] = u(x), 

where r is an arbitrary stopping time with respect to U. Hence, V(x) < u{x) for 
x > a. In particular, if x > b then x < V(x) < u{x) — x and so u{x) = V(x) = 
Ex[U Ta /\T b ] when x > b. 

For the case when a < x < b, note that from (|2.5[) we get for t > that 

Y TbM = M TaAThM + u(x) 

and since 

a — J < Y TbAt <b + J, t>0, 
the optional sampling theorem applies again and we obtain u{x) = E z [y T J. Finally, 
the fact that Y Tb = U TaATb gives us u(x) = E x [U TaATb ] < V(x) and the proof is 
complete. □ 



3. Numerical solution of the free boundary value problem 

We have not been able to give a rigorous proof of the existence and uniqueness 
of the solution (u, b) of the free boundary value problem (|1.4p . We therefore resort 
to a numerical solution by means of the finite element method. However, at the end 
of this section we will show that we have strong computational evidence for both 
existence and uniqueness for (|1.4[) . In order to achieve this we first show rigorous 
existence and regularity results for the boundary value problem (|2.1j) and rigorous 
convergence estimates with explicit constants for the finite element approximation. 

3.1. The boundary value problem. We begin by transforming the free bound- 
ary value problem (|1.4[) to a problem with homogeneous boundary values. Set 
v{x) = u{x) — x and use /_ V^iu) dy = to get 

— 7}<t 2 v"(x) + fixv'(x) 

/oo 
(v(x + y)- v(x))ip(y) dy = -fix, x G (a, 6), 

v(x) =0, x £ (a, 6), 
v'(b) = 0. 

Introducing the operators 

Cv(x) = —t;<t 2 v" (x) + fixv'(x), 



/oo 
(y(x + y)- v(x))<p(y) dy, 
-oo 



our approach will be to first solve the boundary value problem 

Cv —Iv = f, x £ (a, b). 
1 ' ; v(x) = 0, x#(a,b), 

with f(x) = —fix, and then for fixed a < find b > a such that v'(b) = 0. 

To solve (|3. 2[) we follow a standard approach based on a weak formulation and 
Frcdholm's alternative. We denote by (•, •) and || • || the standard scalar product and 
norm in £2(0, b), and we denote by H k (a, b) and i?d(a, b) = {v 6 H 1 (a, b) : v(a) = 
v(b) = 0} the standard Sobolcv spaces. We denote the derivative Dv = dv/dx. We 
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choose v i— » \\Dv\\ to be the norm in Hq(cl, b), which is equivalent to the standard 
ii^-norm. We extend functions v 6 L2(a,b) by zero outside (a, 6) so that Iv is 
properly defined. We define bilinear forms 

Ajc(u, v) = / (gcr u'(x)v'(x) + fj,xu'(x)v(x)j dx, u,v G H (a,b), 

J a 

(3.3) f b 

Ax(u,v) = / Iu(x)v(x)dx, u, v & L 2 (a, b) 



A(u, v) = Ajr(u, v) - A x (u, v). 
Since J™ (p(y) dy = 1, (p(—y) = <p(y), and v(x) = for x ^ (a, 6), we also have 



r-oo 
~oo 



(3.4) Xv(x) = A / ip(x - y)v(y) dy - \v(x), v € L 2 (a,b). 

J a 

The convolution operator Iiw(x) = J_ <p(a; — y)v(y) dy is bounded in £2(0, &) 
with constant c = tp(y) dy = 1 by Young's inequality. Hence, 

(3.5) <2A|H|, veL 2 (a,b), 

(3.6) pXwII < 2A||D«||, v€Hfa,b), 
and 



-Ai^u) > A(jM| 2 - > 0, L 2 (a,6). 



Hence, 



|A(u,«)| < i ( r 2 || J D U .||||^||+ M max(|a|,|6|)!|^|||| W || + 2A|| U |||| W | 
< ci||Dtt||||Dv||, u,v G H^{a,b), 
ci = \<j 2 + c 2 (/imax(|a|, + 2Ac 2 ), 
where we also used Poincare's inequality 

(3.7) ||v|| < ca||£>v||, ueftfM), pa = (6-a)/7r. 
By integration by parts we obtain 

4c(«,«) = ±a 2 ||z^|| 2 - i Ai ||«|| 2 , v G ^(0,6), 
so that A(-, •) is bounded and coercive on Hq (a, b): 

(3.8) \A(u,v)\ < dWDuWWDvW, u,v G H^(a,b), 

(3.9) A(w,w) > iff 2 !!^!! 2 - i/i||«|| 2 , uef#(a,&). 
We say that w G Hq(ci, b) is a weak solution of (|3.2|) if 

(3.10) A(«,0) = (/,0) V^eJ^M). 
We also use the adjoint problem: find w <E H^a, b) such that 

(3.11) A(<f>,w) = (cf>,g) V^eftfM)- 
The strong form is (note that I is self-adjoint in L2(a, b)) 

C*w(x) — Iw{x) = g(x), x G (a,b), 
w(x) =0, x $ (a, b), 

where 

C*w(x) = —^a 2 w"(x) — fj,xw'(x) — fj.w(x). 
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We may now prove the existence and uniqueness of a classical solution of ()3.2[) . In 
principle this follows from the general theory in [5] , but we present a self-contained 
proof, with explicit constants, for the simpler situation that we consider here. The 
theorem also provides results necessary for the analysis of the finite element method. 

Theorem 3.1. The boundary value problem (|3.2[) has a unique weak solution v G 
Hq(o, b) for every f G Li2(a,b). The solution belongs to H (a, b) and there is a 
constant C3 such that 

(3-13) \\D 2 v\\ <c 3 \\f\\. 

Moreover, if f(x) — —fix, then the solution is classical, v G C 2 ([a, b}). Similarly, 
the adjoint problem (|3.12|) has a unique weak solution w G Hq(cl, b) for each g G 
Li2(a,b), which belongs to H 2 (a,b) and 

(3-14) \\D 2 w\\ <c 3 \\g\\. 

Proof. The proof is a standard argument as presented, for example, in [4J Ch. 6] 
for elliptic PDEs. The only difference is that that the lowest order term in A(-, •) 
is defined by means of an integral operator, but the crucial properties (|3.8j) . p.9j) 
are the same. 

We first show that weak solutions are regular. We use a regularity result for 
elliptic problems (see [4j p. 323]): If v is a weak solution of 

Lv{x) = g{x), x G (a, b); v(a) = v(b) = 0, 

and if g G H k (a,b) for some k > 0, then v G H k+2 (a,b). A weak solution v G 
Hl(a, b) of (HI2J) satisfies this with g = f +lv, where by (JHHJ) , (JSH) Iv G if 1 (a, b). 
For / G Li2(a,b) we conclude that v G H 2 (a,b). If / G H 1 {a, b), then we have 
v G H 3 (a,b) and by Sobolev's inbedding v G C 2 ([a,b]). In particular, a weak 
solution is classical when f(x) — and f(x) — —fix. Analogous regularity results 
hold for the adjoint problem. 

Now we can prove existence. Let 

A^u, v) = A(u, v) + v). 

By the Lax-Milgram lemma we know that the shifted problem 

A li {u,4>) = {g,4>) V<£ G Hl[a,b), 

has a unique solution u G Hq (a, b) for each g G £2(0, b). This defines the bounded 
linear operator A^ 1 : £ 2 (a, b) — > Ho(a, b) by u = A~ x g. The equation (|3.10p is now 
equivalent to 

or v — Kv = h, where h = A~ x f and where K = ^fiA^ 1 : L^ia, b) — > £2(0, b) is a 
compact operator, because Hy(a,b) is compactly inbedded in L,2(a,b). 

By the Frcdholm alternative we know that the latter equation is uniquely solvable 
for every h G £2(0, b) if and only if the corresponding homogeneous equation has 
no non-trivial solution. But a non-trivial solution of v — Kv = would be a weak 
solution, and hence a classical solution, of (|3.2|) with / = 0. 

Then we can apply the maximum principle for classical solutions of (|3.2I) , see [5j 
Theorem 3.1.3]. It says that if a classical function satisfies (C —I)u < in (a, b), 
then max[ a £,] u = maxiy a b j u. (The maximum principle for the integro-differential 
equation is proved in the same way as for the differential equation after noting that 
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— Xu(xo) > if u has a maximum at xq.) We conclude that that the homogeneous 
equation has no non-trivial solution and therefore ()3.2|) has a unique weak solution 
for every / £ L2(a,b). By the Fredholm theory the adjoint problem (|3.12l) is then 
also uniquely solvable for all g S L2 

Finally, we prove the bounds (|3.13|) and (|3.14[) . Let u = A^ 1 f and iu = (A*)~ 1 g 
denote the solution operators of (|3.2|) and ()3.12j) . respectively. 

Let / € Hq(ci, b). Then v = A^ 1 f is classical and the maximum principle gives 

( 3 - 15 ) ll u IU«(o,6) < C4||/||l, 00 ( a ,6). 

In order to compute the explicit constant we briefly recall the proof. Let 

( e l(b-a) _ e7 (x-o) x < ^ 

10, .t > 6, 

where 7 > is chosen so that that A(f> > 1 in (a, b). Then u(x) = H/Hz^a, b)4>{ x ) 
satisfies Au > \\f\\L^(a,b) > / = Av in (a,b) and u > = v outside (a,b), so 
that the maximum principle gives maxuM^ — u) = maxjy ( 0) f,) (u — it) = 0, that 
is, u > v in [a, 6]. Hence u < ||0|U oo (a,&)||/|U oo (a,&) in Tnc lower bound 

v > -||^||i 00 (o,6)||/||i 00 (a,6) is obtained in a similar way and so we get 

M\l M ) < 11011^,6)11/11^,6) < e^^WfU^y 
To determine 7, let x € (a, b) and compute 

pb—x 

-14>{x) = \e~< {x - a) / (e™ - l)cp(y) dy 



6 — X" 



Hence, 



A4>(x) > (ia 2 7 2 - fjbj - X)e^ x - a ^ > 1, z e (a, 6), 



if ^cr 2 7 2 — /x&7 — A > 1, that is, if 



V cH 

Then we conclude that (|3.15[) holds with C4 = e^ b_a ) . 

Hence, since ||u|| < (&-a)*||u|Uoo(a,6) and ||/|| Loo (a,b) < (6— a)i , we obtain 
the bound 

H = ll^/ll < c 5 \\Df\\ V/effjJM), c 5 = (b-a)c 4 . 
By duality we conclude 

\U- 1 r\\B(L 2 ,H-i) = U- 1 \\B ( H^L 2) <C 5 - 

Hence 

(3.16) \\wWh-! = IK^*)-^!^-! = IK^- 1 )^!^-! < c 5 |M| V.9 G L a (o,6), 
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where H- 1 {a,b) = (H^(a,b))* and 

\ w\ h-i = sup jjr-rrr. 

Recall that v H> \\Dv\\ is the chosen norm in Hq (a, b). By using <f> = w € H^(a, b) 
here we obtain 

(3.17) \\w\\ 2 <\\w\\ H -4Dw\\. 

We take <j> = w in the adjoint equation (|3.1ip and use coercivity (|3.9p . the inequality 
2afo < ea 2 + e^fr 2 , and (|3~T7|) to get 

±o- 2 j|L>w|| 2 < A(w,w) + ^AH\ 2 < IMilMI + i/i||w'i| 2 

< i//- 1 y| 2 + /illHI 2 < ^'hW 2 + v\H\h-4Dw\\ 

<^-^ 9 \\ 2 + ^a-^w\\ 2 H ^ + l<r 2 \\Dw\\ 2 . 
With (f3~T6l> this leads to 

\\Dw\\ 2 < 2<j- 2 ^ 1 \\g\\ 2 + 4cr-V 2 |HiH-i 
<(2 < 7-V 1 +4a-V- 2 C Dll.9|| 2 
and with Poincare's inequality (|3.7p . 

||w|| < cb||IM| < c 2 (2cr~ 2 At _1 +4cr-V _2 c|)i||ff||. 

Hence 

IIMT^IHIIHI <C6NI V. 9 ei 2 (a,6), 

/o io\ 

c 6 = c 2 (2(T~ 2 n~ 1 + 4<j- 4 n- 2 cl)i. 
By duality in L 2 we also have 

(3-19) IMI = M~V||<cb||/|| V/Gi 2 (a,6). 

In order to bound Z? 2 u we recall that v G H 2 (a, b). Hence it satisfies (|3.2|) strongly, 
so that with (|3.5[) wc obtain 

\a 2 \\D 2 v\\< i x\\xDv\\ + \\Iv\\ + ||/|| 

< M max(|a|,|6|)||Di;|| +2A||v|| + ||/|| 

< ^maxdaUftDHl^wll'lHI* + 2A||t>|| + ||/|| 

< ±a 2 ||i} 2 V || + (2A + a- V max(|a|, |6|) 2 )||«|| + ||/||. 

Hence, 

\\D 2 v\\ <c 7 \\f\\+c % \\v\\, 

c 7 = 4cr~ 2 , c 8 = 4cr~ 2 (2A + i^i + a^ 2 ^ 2 max(|a|, |6|) 2 ). 

In the last step we replaced 2 A by 2 A + fi in c$, so that the same result holds also 
for the adjoint equation (I3.12[) . Using also (|3.18[) and (|3.19[) we finally conclude 

\\D 2 v\\<c 3 \\f\\, \\D 2 w\\<c 3 \\g\\, 
c 3 = c- + c 6 c 8 . 

This completes the proof. □ 
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3.2. The finite element method. The finite clement method is based on a family 
of subdivisions Tk of the interval [a, b] parametrized by the maximal mesh size h. 
Each mesh is of the form 

Th '■ a = xq < x\ < ■ ■ ■ < Xj-i < xj < ■ ■ ■ < x _/y = b, h = max (xj — Xj-i). 

j=l,...,N 

We introduce the space Vh C H${a, b) consisting of all continuous functions that 
reduce to piecewise polynomials of degree < 1 with respect to Th- See Ch. 5] or (TJ 
Ch. 1]. Then there is an interpolator Ih ■ C([a, b]) —> Vh such that Ihu{xj) = u(xj), 
j = 1, . . . , N, and 

(3.20) \\D(u-I h u)\\ Lp{aM < h* + $\\D 2 u\\, ueH 2 (a,b)nH^(a,b), p = 2,oo. 
To prove this we use the identity 

D(u - Ihu)(x) = hj 1 (y! \x) - u '{%))) dy = hj 1 / u"(z)dzdy, 

for x G (xj—i,Xj) and with hj = Xj — Xj-i, which yields 

\D(u - I h u)(x)\ <h]\\D 2 u\\ L2ix3 _ uXj) < hi\\D 2 u\\, x G (xj-t^j). 
This proves the case p = oo and for p = 2 we have 

N 

\\D(u I h u)\\ 2 < ^II^IIW^) < h 2 \\D 2 u\\ 2 . 
j=i 

The finite element problem is based on the weak formulation in (|3.10[) : find 
Vh £ Vh such that 

(3.21) A(vH,<i> h ) = (J,<t> h ) V(f>h€V h , 

where A(-,-) is defined in (|3.3[) with the integral operator computed as in p.4p . In 
the following theorem we prove convergence estimates with explicit constants. 

Theorem 3.2. Let v be the solution of Q3.2[) as in Theorem \3.1\ There is ho = 
o~/(2 2 /is C1C3) such that, for h < ho, p.2ip has a unique solution Vh € Vh and 

(3.22) \\v-v h \\<4c 2 c 2 a- 2 h 2 \\f\\, \\D(v - v h )\\ < 4 Cl c 3 o-- 2 h\\f\\. 

Proof. We adapt an argument from [11] , Let e = v — Vh denote the error. By 
subtraction of (|3.21|) and (|3.10|) with = <\>h € Vh C Hq(o, b) we get 

(3.23) A(e,c(> h )=0 V0 ft G V h . 

Consider the adjoint problem p. lip with g = e and solution w = (A*)~ 1 e. With 
4> = e this yields 

||e|| 2 = A(e,w) = A{e,w-I h w) < Cl \\De\\ \\D(w - I h w)\\ 
< Ci||Z)e||/i||£) 2 u;|| < c lC 3/i||L>e||||e||. 

Here we used (pH3)) . (j3~8| . (j3T20j) . and p~14j) . We conclude 

(3.24) ||e|| < c lC3 h\\De\\. 
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In view of (|3.23|) we have A(e, e) = A(e, v — vt) = A(e, v), so that by (|3.9|) and 
(13341) . 

^ 2 H^|| 2 < A(e, e) + i M |N| 2 = ^(e, «) + |mIN| 2 

< ci||De||||Du|| + \iic\clh 2 \\De\\ 2 . 

Hence, for h < h sufficiently small (ft, 2 , = a 2 / (2fic 2 c 2 )), we have 

||-De|| < c 9 ||L>v||, c 9 = 4cicr -2 . 

Now if / = in (|3.I0p and p.2ip . then v = by uniqueness, and hence e = 0, so 
that Vh — 0. This means that we have uniqueness for the finite element problem 
(13.211) . But this is an equation in a finite dimensional space so existence also follows. 
Therefore, (|3.21[) has a unique solution for all / <G £2(0, b) if h < ho- 

In order to prove the error estimate (|3.22[) we return to (|3.25[) but use A(e, e) = 
A(e, v — = A(e, v — Ihv) instead: 

\(J 2 \\De\\ 2 < A(e,e) + ^i\\e\\ 2 = A(e,v- I h v) + i/i||e|| 2 
< Cl \\De\\\\D{v - I h v)\\ + ^c 2 lC lh 2 \\De\\ 2 , 
and conclude, for h < ho, 

||£>e|| < cg\\D(v - I h v)\\, c 9 = Ac 1( t- 2 . 
Hence, by (pT20j) . (j3~T3l) . and (j3~24| . 

\\De\\ < c 9 h\\D 2 v\\ < c 9 c 3 h\\f\\ = Ac lC3 a- 2 h\\fl 
\\e\\<c lC:i h\\De\\<4c 2 c 2 a- 2 h 2 \\f\\, 
which is (pH!2"]) . □ 
We finish by proving the pointwise convergence of the derivative. 

Corollary 3.3. Assume that each finite element mesh Th is uniform, that is, Xj — 
Xj-i = h for j = 1, . . . , N. Then, for h < ho as in Theorem \3.2l we have 

\v'(b)-v' h (b)\ <c 10 h*\\f\\, c 10 = 2 + 4 Cl c 3 o-- 2 . 
Proof. We use the inverse inequality 

||£%IUooO,i>) ^ h-i\\D<t> h \\, <j> h G V h . 
To prove this we note that 



D(j) h (x) = h I D(j) h (y)dy, x <G Xj), h = x 3 ■- Xj-i, 

which yields 

\D<j) h (x)\ < h-i\\D<t> h \\ L2{Xj _ uXj) < h-^\\Dtj> h \\, x e (xj-^xj). 
Hence, by (j330|) and (j3~22|) . 

ll-De|| ioo ( 0i6 ) < \\D(v -I h v)\\ L „( a ,b) + \\D(I h v - v h )\\ Lx{a , b) 

< \\D(v - I h v)\\ LaoM +h-*\\D(I h v-v h )\\ 

< \\D(v-I h v)\\ LaoM + h-i\\D(I h v-v)\\+h-i\\D(v-v h )\\ 

< 2hi\\D 2 v\\ +h-i\\D(v-v h )\\ < (2 + 4 Cl c 3 0^ 11/11- 
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Therefore 

\v'(b)-v' h (b)\ < (2 + 4 Cl c 3 a~ 2 )hi\\f\\. 

□ 

In particular, with f(x) = —fix, Corollarv l3.3l gives 



b 3 



(3.26) \v'(b) - v' h (b)\ < cnh\ c n = cioMy 

Given numerical values for the parameters a, b, a, /i, A we may now compute 
numerical values for ho and cn. Alternatively, we may conclude that there are 
uniform bounds > h^, c\\ < C\\ for b G [^1,^2] and with the other parameters 
fixed. 

3.3. The free boundary value problem. We use uniform meshes Th with 

b — a 

x j - x j-i = h = — a^: J = 1,...,N. 

Since we want to vary b, we parametrize by N instead of h. Let f(x) = —fix, fix 
a < and let v, vn denote the solutions of (|3.10|) and (|3.21|) for b > a. Define the 
functions 

F(b) = v'(b), F N (b)=v' N (b). 
From p.26|) . we get for a < b\ < 62 

\\F - F N \\ Loo(bub2) < c 12 N~i , N>N , 

(3-27) „ , b 2 -a 

C12 =cii(0 2 -a) 2 , A'o = — ■ 

By writing down the matrix equation for solving the finite element problem (|3.2f I) , 
it is easy to see that, for fixed TV, the function b i-> F^(b) is continuous on (a, 00). 
From (|3.27[) we conclude that b H> F(b) is also continuous on (a, 00). Moreover, 
by a direct consequence of the strong maximum principle and the Hopf boundary 
point principle for our equation (see [H Theorem 3.1.4-3.1.5]), we get the following: 

Lemma 3.4. If a < b < 0, then F(b) < 0. In particular, if (u, b) is a solution to 
the free boundary problem (jl.4p . then b > 0. 



We believe that there exists a unique b > such that F(b) = 0. We are not 
able to provide a rigorous proof of this, but numerical simulations present strong 
evidence in the following way. Assign numerical values to the parameters a, a, \i, A 
and fix a jump density ip. In all our computations, we took tp to be the truncated 
normal distribution with mean zero, variance 7 > and support [—J, J], i.e. 

f(v) = < 7V ^F(2$(j/ 7 )-i) if - J<y<J > 



otherwise, 



where 



= — = / e 2 dy, x g 
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From computations of the boundary value problem (|3.21j) (see Figures l3~Tl and l3~2|) . 
we can find < bi < b 2 and N > Nq such that 

Fjfih) < ~\, F^{b 2 ) > 1 and c 12 N-? < i. 

(The 1/2 and 1/4 may vary if we change the parameters.) From (|3.27[) . we can 
then conclude that 

F(h) < 0, F(b 2 ) > 0, 

Fjv(6i) < 0, F N (b 2 ) > for all N > N. 

Hence, there exists b G (b\, b 2 ) such that F(b) — and for each N > N there exists 
&at G {bi,b 2 ) such that F N (bN) = 0. Moreover, (|3.27p gives us that 

lim F(b N ) = 0. 

Of course, we cannot conclude that b is unique and &at — >■ & as N — > oo. However, 
Figure [3~T1 suggests that 6 is unique and from computations with increasing N, it 
seems like b^ converges, see Table [3~T1 

Wc now discuss whether the properties a) and b) in the statement of Thcorcm ll.il 
hold for a solution (u,b) of (|1.4|) . We have no rigorous proof, but computational 
evidence. The properties a) and b) boil down to 

,6 

(3.28) A / v{y)ip{y — x) dy < fix, for x > b, 

J a 

and v > respectively, where (v, b) solves (|3.1|) . We believe that u > holds for all 
values of the parameters, but computations suggests that (|3.28|) may fail for certain 
parameter values, typically when a is small and A is three or four times larger than 
fi. See Figures 1331 and I3~4l where we check (|3.28[) for (vN,bjsr) instead of (v,b). 



N 


b N 


2000 


0.0572939 


4000 


0.0572743 


6000 


0.0572678 


8000 


0.0572653 



Table 3.1. a = -0.1, A = 10, a = 0.2, fi = ^ 5E , 7 = 0.02 and 
J = 0.05. 
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Figure 3.2. The solution (v N ,b N ) when a 
0.2, n = 7 = 0.02 and J = 0.05. 



= -0.1, A = 10, a = 
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FIGURE 3.3. A simulation of (|3~28f when a = -0.1, A = 30, a = 
0.2, n = g^g, 7 = 0.02 and J = 0.05. The condition fails. 
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FIGURE 3.4. A simulation of (|3.28J) when a = -0.1, A = 10, a = 
0.2, /i = g^g, 7 = 0.02 and J = 0.05. The condition holds. 



